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Abstract 

Recent scanning tunneling experiments on Ca2-x^0,xCuO2Cl2 by Hanaguri et al.- observe field- 
dependent quasiparticle interference effects which are sensitive to the sign of the d-wave order 
parameter. Their analysis of spatial fluctuations in the local density of states shows that there 
is a selective enhancement of quasiparticle scattering events that preserve the gap sign, and a 
selective depression of the quasiparticle scattering events that reverse the gap sign. We introduce 
a model which accounts for this phenomenon as a consequence of vortex pinning to impurities. 
Each pinned vortex embeds several impurities in its core. The observations of recent experiments 
can be accounted for by assuming that the scattering potentials of the impurities inside the vortex 
cores acquire an additional resonant or Andreev scattering component, both of which induce gap 
sign preserving scattering events. 

PACS numbers: 74.20.-z; 74.25. Jb; 74.72.-h 
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I. INTRODUCTION 



Fundamental studies of unconventional superconductors are currently hindered by the 
scarcity of direct methods to determine the structure of the superconducting order parame- 
ter. Apart from Josephson junction experiments, few spectroscopic probes provide the valu- 
able information about the phase of the order parameter. In this work, we discuss how phase 
sensitive coherence effects can be studied using scanning tunneling spectroscopy/microscopy 



The key idea is that the evolution of the phase of the order parameter in momentum 
space can be determined from the Fourier transformed fluctuations in the tunneling density 
of states. The sensitivity of these fluctuations to the scattering rates of superconducting 
quasiparticles manifests itself through coherence factor effects. Quasiparticles in a super- 
conductor are a coherent superposition of excitations of electrons and holes. Coherence 
factors characterize how the scattering rate of a superconducting quasiparticle off a given 
scatterer differs from the scattering rate of a bare electron off the same scatterer-^i. Coher- 
ence factors are determined by combinations of the Bogoliubov coefficients and v^, which 
give proportions of the particle and hole components that constitute a superconducting 
quasiparticle. 



The momentum-dependent order parameter = | Ak|e*'^'^'*^ has the same sign as the Bo- 
goliubov coefficient v^, so that studies of scattering rates of quasiparticles with different 
momenta can delineate how the phase of the order parameter 0(k) changes in momentum 
space. 

In studies of unconventional superconductors with spatially varying order parameter, 
scanning tunneling spectroscopy provides a spectroscopic probe with a real space resolution 
at the atomic level. In the past, observation of phase sensitive coherence effects with STM 
has been thwarted by the problem of controlling the scatterers^. An ingenious solution of 
this problem has been found in the application of a magnetic field, which introduces vortices 
as controllable scatterers in a given system^. 

In this work, we develop a framework observation of coherence factor effects with Fourier 



(STS/STM). 



(1) 
(2) 
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Transform Scanning Tunneling Spectroscopy (FT-STS). Using this framework, we analyze 
the recent observations of the coherence factor effects in a magnetic field to develop a phe- 
nomenological model of quasiparticle scattering in a disordered vortex array. 



II. COHERENCE FACTORS IN STM MEASUREMENT 

Scanning tunneling spectroscopy, which involves tunneling of single electrons between 
a scanning tip and a superconducting sample, offers an opportunity to examine how the 
spectrum of superconducting quasiparticles responds to disorder. We now discuss how we 
can extract phase-sensitive information from STM data. 

A. LDOS correlators R"^^^^ and have well-defined coherence factors 

We describe the electron field inside a superconductor by a Balian-Werthammer spinor- 

( ^T(r,r)\ 

^|(r,r) 

where r denotes real space coordinates and r is imaginary time. The Nambu Green's function 
is defined as the ordered average 



^fr, t) 



G„^(r',r;r) = -(T.vl>„(r', r)vl/t (r, 0)). 
Tunneling measurements determine local density of states, which is given by 

p(r, a;) = — Ira Tr ^ [^(r, v]u — i5)] , 

TT 2 



(3) 



(4) 



where G(r', r; z) is the analytic continuation G{y' , r; iuOn) G{r', r; z) of the Matsubara 
Green's function 

G{r',r;tu;n)= G{v\v-T)e'^-^ dr, (5) 
Jo 

with ojn = {2n + 1)'kT. The appearance of the combination ^-y^ in (jll) projects out the 
normal component of the Nambu Green's function 



Tr G(r', r; r) = - J^i^rM^'. r)^Ptir, 0)). 



(6) 



The mixture of the unit and the T3 matrices in this expression prevents the local density of 
states from developing a well-defined coherence factor. We now show that the components 
of the local density of states that have been symmetrized or antisymmetrized in the bias 
voltage have a well-defined coherence factor. The key result here is that 



p(r, uj) ± p(r, —u) = — Im Tr 

71 
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(7) 



In particular, this implies that the antisymmetrized density of states has the same coherence 
factor as the charge density operator T3. 

To show these results, we introduce the "conjugation matrix" C = (T2^2, whose action on 
the Nambu spinor is to conjugate the fields. 



(8) 



effectively taking the Hermitian conjugate of each component of the Nambu spinor. This 
also implies that \1/^C = ^P^. Here Tj are Pauli matrices acting in particle-hole space, for 
example. 



and o"j are Pauli matrices acting in spin space. 





Using ([8]), it follows that 

[CG{T',r;T)CU = -(T,CM/(r',r)M/t(r,0)C)«^ 

= (T.vl/^(r,0)vl/l(r',r)) 
= -G/3a(r,r', -r), 

or, in the matrix notation, 

CG{r,r';T)C = -G^(r', r; -r), 
which in turn implies for the Matsubara Green's function ([5]) 

CG{r,r']iuJn)C = -G^{r' ,r; -iun)- 



(9) 



(10) 



(11) 



For the advanced Green's function, which is related to the Matsubara Green's function via 
analytic continuation, G{r,r' ,iun) G{T,r',z), we obtain 



CG{r,r';u-iS)C = -G^ {r' ,r; -u + iS). 



(12) 



Using this result and the commutation relations of Pauli matrices, we obtain 

1 1 + T3 

p(r, —u) = Im Tr — - — G{r, r; —u + i6) = 

71 2 

= - Im Ti G G^(r, r;uj-i6) G = 

TT 2 

1 1 - Ta 
= — Im Tr ^(r, r;uj — i5). 

TT 2 



(13) 



Finally, we obtain 



p(r, uj) ± p(r, —uj) = — Im Tr 



~|~ G'(r, r; (u; — i(5) ± — — ^ G'(r, r; u; — i5) 



Im Tr 



vr 



G(r, r; cij — i6) 



T3 



(14) 



B. Coherence factors in a BCS superconductor, T-matrix approximation 

Next, applying this result to a BCS superconductor, we show that in the t-matrix ap- 
proximation the coherence factors that arise in the conductance ratio ^(q, V) are given by 
the product of the coherence factors associated with the charge operator and the scattering 
potential. 

T-matrix approximation^"^ allows to compute the Green's function in the presence of 
multiple scattering off impurities. In terms of the bare Green's function G(k, uj) and the 
impurity t-matrix t(k, k'), the full Green's function is given by 

G(k, k', uj) = G(k, uj) + G{k, uj)i{k, k')G{k', uj) = G{k, uj) + 6G{k, k', uj). (15) 



Using this expression, we obtain for the Fourier transformed odd fluctuations in the tunneling 
density of states 



Sp°'^'^{q,u) = r^lm I Tr 



1 

2^" 



T^3'^G'k+,k_ - iS) 



271 



-Im / Tr 



TsGk- (t^ - iS)i{q, k)G'k+ {uj - i6) 



(16) 
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The Fourier transformed even fluctuations in the tunnehng density of states 



iq, ^} = J Tr 



2tt 



-Im / Tr 



(17) 



For scattering off a single impurity with a scattering potential ?7(k, k'), the t-matrix 
t(k, k') denotes the infinite sum 



t(k, k') = U{k, k') + J2 U{K k")G{k", u)U(k", k') + ... = 

k" 

= f/(k, k') + J2 U{K k")G{k", uj)i{k", k') 



(18) 



Working in the Born approximation, which is equivalent to taking only the first term 
in the series (1181) . we derive the expressions for the coherence factors associated with some 
common scattering processes that arise in the even and odd density-density correlators 
R^^'^'^{q,V) and -R°'^'^(q, y) in a BCS superconductor (see Table 1). We use the following 
expression for the BCS Green's function for an electron with a normal state dispersion ek 
and a gap function Ak: 



Gk{uj) = [uj- CkTs - AkTi] \ 



(19) 



t(q, k) is the scattering t-matrix of the impurity potential, and k± = k ± q/2. If the 
scattering potential has the t-matrix given by t(q, k) = T3(q) t^, corresponding to a weak 
scalar (charge) scatterer, the change in the odd part of the Fourier transformed tunneling 
density of states becomes S p°^i^^{q, uj) = T^iq) A°f;^;^^(q, ^) with 



27r .Ik 



+ ek+ek_ - Ak, Ak_ 



EL 



(20) 



where i?k = [Ck + Ak]^ is the quasiparticle energy. Expressed in terms of the Bogoliubov 
coefficients and fk, given by uKvl) = |(1 ± tk/Ek), the expression under the integral in 
(!20l) is proportional to (m+u_ — t>+t>„)^. 

Fluctuations in the even part of the Fourier transformed tunneling density of states due 
to scattering off a scalar impurity are substantially smaller, i?^^^;^j,(q, u;) <^ R°f^i^j.{q, u) , 



where R 



,even(pdd) 
scalar 



q,^) is defined by V^cX.lq,^) = ^3(q) KZLH^) with 



Expressed in terms of the Bogoliubov coefficients Wk and f k, the expression under the integral 
in fl2Tl) is proportional to + — and is, therefore, small for the 

nodal quasiparticles involved, |A^^^^^(q, cu) | ^ |A^^^;^^(q, ct;)|. Thus, scattering off a weak 
scalar impurity contributes predominantly to odd-parity fluctuations in the density of states, 
^odd(q^ 1/). 

In a second example, consider scattering off a pair-breaking "Andreev" scatterer with the 
t-matrix given by t(q, k) = Ti(q, k) ri. Here the change in the even and odd parts of the 
Fourier transformed tunneling density of states are (5p^'^"^°'^'^^(q, cj) = A'^J^^"'^°'^'^\q_, uj) with 



1 



27r 
1 



k 



2;(Ak+ + A 



Ar"(q,^) = 7^Im / ri(q,k) — ::^-|±^— (22) 



Ck, Ak_ + ek_A 



Z=UJ — t< 
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In terms of the Bogoliubov coefficients Wk and fk, the expressions in square brackets in 
A'^j^^"' {q, u) and A'^'^{q,uj) are proportional to + i;+t>_)(u+t>_ + and (u^u_ — 

f (-u+f _+f respectively. For the nodal quasiparticles involved, the latter expression 
is substantially smaller than the former, |A^°'°'(q, u;)| ^ |A^^"(q, cj)!. Thus, scattering off 
an Andreev scatterer gives rise to mainly even parity fluctuations in the density of states. 

We summarize the coherence factors arising in i?'^''^"(q, V) and i?°'^'^(q, V) for some com- 
mon scatterers in Table 1. The dominant contribution for a particular type of scatterer is 
given in bold. 

TABLE I: Coherence factors C(q) in i?^^'^"(q, and R°'^'^{(\,V) for some common scatterers. 



T-matrix 


Scatterer C(q) in i?^^""(q, V) 


C(q) in R°'^'^{(i, V) 


Enhanced 


qi Enhances 


T3 


Weak Scalar [uu' + vv'){uu' — vv') 


(uu' - vv')2 


2,3,6,7 


No 


fj • in 


Weak Magnetic 





None 


No 


i sgn io 1 


Resonant (uu' + vv')-^ 


{uu' + vv'){uu' — vv') 


1,4,5 


Yes 


Tl 


Andreev (uu' + vv')(uv' + vu'^ 


{uu' — vv'){uv' + vu') 


1,4,5 


Yes 



From Table 1, we see that the odd correlator R°'^'^{(\,V) is determined by a product of 
coherence factors associated with the charge operator and the scattering potential, while the 
even correlator i?^^*^"(q, V) is determined by a product of the coherence factors associated 
with the unit operator and the scattering potential. 
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C. Conductance ratio - measure of LDOS 



An STM experiment measures the differential tunneling conductance J^(r, V) at a loca- 



dl 

iiiuuuifciiiut; 

tion r and voltage In a simplified model of the tunneling 



dl r 

— (r,y)oc J ^diu[-f{LU-eV)] J dr,dr2M{r,,r)M''{r2,r)A{r2,r,,u;), 



(24) 



where A{r2,Ti,u) = \lm G{y2,Ti,uj — i5) is the single electron spectral function and f{uj) 
is the Fermi function. Here ri, T2 and r are the two-dimensional coordinates of the incoming 
and outgoing electrons, and the position of the tip, respectively. M(ri, r) is the spatially de- 
pendent tunneling matrix element, which includes contributions of the sample wave function 
around the tip. 

Assuming that the tunneling matrix element is local, we write M(ri, r) = M(r)5*^^)(ri— r), 
where M(r) is a smooth function of position r. In the low-temperature limit, when T — > 0, 
the derivative of the Fermi function is replaced by a delta-function, —fioj — eV) = 6{uj — eV). 
With these simplifications, we obtain 

^(r,\/)cx|M(r)|V(r,l-), (25) 

where p(r, V) = A{r, r, V) is the single-particle density of states. In the WKB approach 
the tunneling matrix element is given by |M(r)p = e~^^^^^ with 7(r) = f^^^^ dx- /^HJ^^tM. 



s(r) 



^ ■ -y/2m^^(r) , where s(r) is the barrier width (tip-sample separation), i/j{r) is the barrier 
height, which is a mixture of the work functions of the tip and the sample, m is the electron 
mass^i^. Thus, the tunneling conductance is a measure of the thermally smeared local density 
of states (LDOS) of the sample at the position of the tip. 

To filter out the spatial variations in the tunneling matrix elements M(r), originating from 
local variations in the barrier height and the tip-sample separation s, the conductance ratio 
is taken: 

7(r V) = ^^"-+^^ = +V) ^ Poi+V) + Spir,+V) 
^' ' §{r,-V) p{r,-V) po{-V) + 6p{r,-Vy ^ ^ 

For small fluctuations of the local density of states, Sp{r, ±V) <^ po{±V), Z{r, V) is given 

by a linear combination of positive and negative energy components of the tunneling density 

of states, 

Z(..t/)==Z„M[l + M!i±^_M^l (27) 

L Po{+V) Po{-V) J 
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with Zo{V) = ^Ji^v) ■ '^^^ Fourier transform of this quantity contains a single delta function 
term at q = plus a diffuse background, 

Z(q, V) = MV)(2.m<t) + Z„(V) _ M^l . (28) 

Interference patterns produced by quasiparticle scattering off impurities are observed in the 
diffuse background described by the second term. 

Clearly, linear response theory is only valid when the fluctuations in the local density of 



states are small compared with its average value, 5p{r, ±Vy <^ po{±VY. In the clean hmit, 
this condition is satisfied at finite and sufficiently large bias voltages |V| > 0. At zero bias 
voltage V ^ 0, however, the fluctuations in the local density of states become larger than 
the vanishing density of states in the clean limit, \dp{r, ±V)\ > po{±V), and linear response 
theory can no longer be applied. 

At flnite bias voltages, \V\ > 0, fluctuations in the conductance ratio Z{q,V) are given 
by a sum of two terms, even and odd in the bias voltage: 

Z(q, y)|q^o = ZoiV) hp^^^^iq, V)i—^ - — + Sp'^'^id, V)i—^ + ^ 



Po{+V) po{-V)' ^ ''po{+V) po{-V) 

(29) 

where 5pe^e"H'')(q, V) = (5p(q, +V) ± 5p(q, -V))/2. 

Depending on the particle-hole symmetry properties of the sample-averaged tunneling 
density of states Po{V), one of these terms can dominate. For example, if at the bias voltages 
used, the sample-averaged tunneling density of states Pq{V) is approximately particle- hole 
symmetric, po{—V) ^ po{+V) — po{V), then Z(q,V) is dominated by the part of LDOS 
fluctuations that is odd in the bias voltage V, 

Z(q, V)\^^o ^ Zo{V) ^Sp'^^^id, V). (30) 

In general, when we average over the impurity positions, the Fourier transformed fluc- 
tuations in the tunneling density of states, 5p(q, V), vanish. However, the variance in the 
density of states fluctuations is non-zero and is given by the correlator 



R{q,V) - 5p{q,V)6p*{-<i,V). (31) 

Defining 



we obtain that for q 7^ 



4|^o(K)P 

Pliv) 



(33) 



D. Observation of coherence factor effects in QPI: coherence factors and the octet 
model 

In high-Tc cuprates the quasiparticle interference (QPI) patterns, observed in the Fourier 
transformed tunneling conductance dl{q, V)/dV oc p(g, V^), are dominated by a small set of 
wavevectors ^1-7, connecting the ends of the banana-shaped constant energy contours^i^i^. 
This observation has been explained by the so-called "octet" model, which suggests that 
the interference patterns are produced by elastic scattering off random disorder between the 
regions of the Brillouin zone with the largest density of states, so that the scattering between 
the ends of the banana-shaped constant energy contours, where the joint density of states is 
sharply peaked, gives the dominant contribution to the quasiparticle interference patterns. 

In essence, the octet model assumes that the fluctuations in the Fourier transformed 
tunneling density of states are given by the following convolution: 



While this assumption allows for a qualitative description, it is technically incorrect^ii^,, for 
the correct expression for change in the density of states involves the imaginary part of a 
product of Green's functions, rather than a product of the imaginary parts of the Green's 
function, as written above. In this section, we show that the fluctuations in the conductance 
ratio at wavevector q, given by Z(q, V^), are, nevertheless, related to the joint density of 
states via a Kramers-Kronig transformation, so that the spectra of the conductance ratio 
Z(q, V) can still be analyzed using the octet model. 

As we have discussed, fluctuations in the density of states 5p(q, V) are determined by 
scattering off impurity potentials and have the basic form ( JT6l) . This quantity involves the 
imaginary part of a product of two Green's functions, and as it stands, it is not proportional 
to the joint density of states. However, we can relate the two quantities by a Kramers-Kronig 
transformation, as we now show. 

We write the Green's function as 





(34) 



10 



where G'^Sui — iS) = jr{Gi^{uj — iS) — Gkl^^ + iS)). Substituting this form in (fT6l) . we obtain 
Sp"''{q,E)= ^ [ Tr[r3 [ dE' [^^J2^'UE)iMG'UE') - [E ^ E']]]{35) 

Jk J 

As we introduce the joint density of states, 

J(q, E, E') = ^ Tr[r,G'UE)i{c^, k)^;^^^')], (36) 

( |35|) becomes 

5p°^\c^, u:)=\jdE' [ J(q, E,E') + J(q, E\E)]. (37) 

The Fourier transformed conductance ratio Z{q, E) given by fl30|) now becomes (for q 7^ 0) 
Z(q, i?) = / rfi?' ^[^(q, E, E') + J(q, E', E)]. (38) 
Substituting the expression for the BCS Green's function ( fT9l) in ( l36l) . we obtain 



J(q, ^, = i e;:^^^[^3(^ + ek-Ta + Ak-ri)t(q, k)(E' + Ck+rg + Ak+Ti)] • 

• [6{E - Ek_) - S{E + ^k-)] [S{E' - Ek+) - 5(^' + ^k+)] ■ sgnE ■ sgnE'{39) 



where i?k± = y ek± + ^k±- Provided both the energies are positive, E, E' > 0, we obtain 

J(q,E,E') = Ek,k. C(ki,k2)5(E - E^,)5{E' - i?kJ5(')(ki - k2 - q), (40) 

where the coherence factor is 

C(ki, ka) = ^— ^Tr[r3(^ + e^,n + Ak,ri)£(ki, k2)(^' + ek.rg + Ak.ri)]. (41) 
4 -fcki-c^ka 

Now the fluctuations in the conductance ratio at wavevector q are given by: 

/jTpi r 
^ j dkidk2C(ki, k2)5(E - E^,)5{E' - (3E^,)5^^\\^, - k2 - q).(42) 

Thus, the fluctuations in the conductance ratio 2'(q, E) are determined by a Kramers-Kronig 
transform of the joint density of states with a well-defined coherence factor. 

Conventionally, coherence factors appear in dissipative responses, such as fHU]) . The 
appearance of a Kramers-Kronig transform reflects the fact that tunneling conductance is 
determined by the non-dissipative component of the scattering. The validity of the octet 
model depends on the presence of sharp peaks in the joint density of states. We now argue 
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that if the joint density of states contains sharp peaks at well-defined points in momentum 
space, then these peaks survive through the Kramers-Kronig procedure, so that they still 
appear in the conductance ratio 2'(q, E) with a non-Lorentzian profile, but precisely the 
same coherence factors. We can illustrate this point both numerically and analytically. Fig. 
1 contrasts joint density of states with the Fourier transformed conductance ratio 2'(q, E) 
for scattering off a weak scalar impurity, showing the appearance of the "octet" scattering 
wavevectors in both plots. Similar comparisons have been made by earlier authors^^ii^. 

Let us now repeat this analysis analytically. Suppose J(q, Ei, E2) (HOj) has a sharp peak at 
an octet q vector, q = qj (i = 1 — 7), defined by the delta function J(q, Ei = E , E2 = E) = 
Cj(5''^^(q — qj), where Ci is the energy-dependent coherence factor for the ith octet scattering 
process. When we vary the energy E2 away from E, the position of the characteristic octet 
vector will drift according to 

q,(Ei, E2) = q^{E) - VeMEi - E) + V eME2 - E), (43) 

where V^i^q^ = ^ni(i) and V£;2qi = ^n2(«) are directed along the initial and final quasi- 
particle velocities, and f a is the quasiparticle group velocity. Carrying out the integral over 
E' in fH2l) we now obtain 



Z(q, E) ^ [ dE' kq - q.(E) -"^{E'- E)) + <5(q - c^,{E) + ^{E' - E)) 

J E - E' I va va 

jz—TT-^ ((q - - . ^v/ ((q - q^)^^)] , (44) 



= Q 

where 

(q- qi)iii,2 = (q- q*) ■ ni,2(i) 

denotes the component of (q — q^) parallel to the initial/final quasiparticle velocity and 

(q - qi)±i,2 = (q - q*) ■ [z x ni,2(0]- 

denotes the component of (q — q^) perpendicular to the initial/final quasiparticle velocity, 
where z is the normal to the plane. Thus, a single sharp peak in the joint density of states 
produces an enhanced dipolar distribution in the conductance ratio Z{ci,E), with the axes 
of the dipoles aligned along the directions of the initial and final quasiparticle velocities. 
The above analysis can be further refined by considering the Lorentzian distribution of the 
quasiparticle interference peaks, with the same qualitative conclusions. 
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(a) (b) 

FIG. 1: (Color online) Observation of coherence factor effects in the squared joint density of states 
I J(q, y, and in the squared Fourier transformed conductance ratio \Z{<\^ Fig- (a) shows 
the squared joint density of states |J(q, 1/)^ at the bias voltage V = Ao/2, Fig. (b) shows the 
squared Fourier transformed conductance ratio |Z(q, 1/)^ produced by a weak scalar scattering 
potential t(q) = fs. Red lines label the positions of the sign-reversing q-vectors q = g2,3,6,7> where 
weak scalar scattering is peaked. Blue lines label the positions of the sign-preserving q-vectors 
q = (J'1,4,5, where weak scalar scattering is minimal. 



To summarize, the conductance ratio 2'(q, E) is a spectral probe for fluctuations in the 
quasiparticle charge density in response to disorder. Z'(q, E) is characterized by the joint 
coherence factors of charge (t^) and the scattering potential. Provided the original joint 
density of states is sharply peaked at the octet vectors q^, i = 1 — 7, the conductance ratio 
Z{q, E) is also peaked at the octet vectors qj, z = 1 — 7. 



III. MODEL FOR QUASIPARTICLE INTERFERENCE IN VORTEX LATTICE 

Next, we discuss the recent experiments by Hanaguri et al.^ on the underdoped cuprate 
superconductor calcium oxychloride, Ca2-xNaxCu02Cl2 (Na-CCOC), which have success- 
fully observed the coherence factor effects with Fourier Transform Scanning Tunneling Spec- 
troscopy (FT-STS) in a magnetic field. The main observations are: 
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• A selective enhancement of sign-preserving; depression of sign-reversing 
scattering events. In a field, Hanaguri et al.- observe a selective enhancement of the 
scattering events between parts of the Brillouin zone with the same gap sign, and a 
selective depression of the scattering events between parts of the Brillouin zone with 
opposite gap signs, so that the sign-preserving q-vectors qi,4,5 are enhanced, and the 
sign-reversing q-vectors q2,3,6,7 are depressed. 

• Large vortex cores with a core size ^ ~ 10a of order ten lattice constants. 
Experimentally, vortex cores are imaged as regions of shallow gap^. The figure 
C, ~ 10a is consistent with magnetization and angular resolved photoemission (ARPES) 
measurements^^ . 

• High momentum transfer scattering involving momentum transfer over a large 
fraction of the Brillouin zone size at ^4^5 ^ kp- A paradoxical feature of the observa- 
tions is the enhancement of high momentum transfer q n/a scattering by objects 
that are of order ten lattice spacings in diameter. The enhanced high momentum 
scattering clearly reflects sub-structure on length scales much smaller than the vortex 
cores. 

• Core-sensitivity. Fourier mask analysis reveals that the scattering outside the vor- 
tex core regions differs qualitatively from scattering inside the vortex core regions. 
In particular, the enhancement of the sign-preserving scattering events is associated 
with the signal inside the "vortex cores" , whereas the depression of the sign-reversing 
scattering events is mainly located outside the vortex regions. 

Recently, T. Pereg-Barnea and M. Franz^^ have proposed an initial interpretation of these 
observations in terms of quasiparticle scattering off vortex cores. Their model explains the 
enhancement of the sign preserving scattering in the magnetic field in terms of scattering off 
vortex cores, provided vortex cores are small with ^ ~ a, as in high temperature superconduc- 
tor Bi2Sr2CaCu208+5 (Bi2212). However, the large vortex core size of Ca2-xNaxCu02Cl2 
is unable to account for the field-driven enhancement in the high momentum scattering. 

Motivated by this observation, we have developed an alternative phenomenological model 
to interprete the high-momentum scattering. In our model, vortices bind to individual 
impurities, incorporating them into their cores and modifying their scattering potentials. 
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This process replaces random potential scattering off the original impurities with gap-sign- 
preserving Andreev reffections off order parameter modulations in the vicinity of the pinned 
vortices. The high-momentum transfer scattering, involved in the selective enhancement 
and suppression, originates from the impurities whose scattering potentials are modified by 
the presence of the vortex lattice. Rather than attempt a detailed microscopic model for 
the pseudo-gap state inside the vortex cores and impurities bound therein, our approach 
attempts to characterize the scattering in terms of phenomenological form factors that can 
be measured and extracted from the data. 

A. Construction of the model 

In the absence of a field, random fluctuations in the tunneling density of states are 
produced by the original impurities. We assume that scattering off the impurities is mutually 
independent permitting us to write the change in density of states as a sum of contributions 
from each impurity 

Spiv, V,B = 0) = J2 SPri^ - V)^ (45) 

j 

where denote the positions of the impurities. If 

rii = original concentration of impurities in the absence of magnetic field, 
then we obtain 

i?(q, V,B = 0)= 5pi{q,V)Sp*i-q,V). (46) 

Next we consider how the quasiparticle scattering changes in the presence of a magnetic 
field. Pinned vortices arising in the magnetic field act as new scatterers. In the experimentii, 
vortices are pinned to the preexisting disorder, so that in the presence of a magnetic field, 
there are essentially three types of scatterers: 

• bare impurities, 

• vortices, 

• vortex-decorated impurities. 
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Vortex-decorated impurities are impurities lying within a coherence length of the center of 
a vortex core. We assume that these three types of scattering centers act as independent 
scatterers, so that the random variations in the tunneling density of states are given by the 
sum of the independent contributions, from each type of scattering center: 

3 I m 

where rj,rj,r^ denote the positions of vortices, decorated impurities and bare impurities, 
respectively. In a magnetic field, the concentration of vortices is given by 

2eB 

nv = concentration of vortices = — — , 

h 

In each vortex core, there will be Ucore = '^j7r(^^/4) impurities, where 7r(^^)/4 is the area of 
a vortex and rii is the original concentration of bare scattering centers in the absence of a 
field. The concentration of vortex-decorated impurities is then given by 

2eB 

udi = concentration of vortex-decorated impurities = Ucorenv = — 7— ^i7r(^/2)^. 
Finally, the residual concentration of "bare" scattering centers is given by 

nj = rii — Udj = concentration of residual "bare" impurities. (48) 
Treating the three types of scatterers as independent, we write 



i?(q, V) = nv 5py(q, V)6p*yi-(i, V) + udi 5pD/(q, V)6p*j,j{-(i, V) 



+ {m-nDi) 6pi{q,V)6p*j{-q,V). (49) 

The first term in accounts for the quasiparticle scattering off the vortices, the second 
term accounts for the quasiparticle scattering off the vortex-decorated impurities and the 
third term accounts for the quasiparticle scattering off the residual bare impurities in the 
presence of the superflow. It follows that 

\Ziq,V,B)\' = ^ \Zyiq,V,B)\' + ^n,,re\Zniiq,V,B)\' + {n,-^n,ore) \Zjiq,V, B)\', 
h h n 

(50) 

where Z{q,V,B) is given by ( l30l) . averaged over the vortex configurations, Zv{q,V, B), 
Z£,j{q,V) and Zj{q,V) are Fourier images of the Friedel oscillations in the tunneling den- 
sity of states induced by the vortices, vortex-decorated impurities and the bare impurities 
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in the presence of the superflow. Our goal here is to model the quasiparticle scattering phe- 
nomenologically, without a recourse to a specific microscopic model of the scattering in the 
vortex interior. To achieve this goal, we introduce Zvi{q, V,B), a joint conductance ratio of 
the vortex-impurity composite, which encompasses the scattering off a vortex core and the 
impurities decorated by the vortex core, 

l^Vll"^ = l^y^ + n'core\ZDl\'^, (51) 

so that we obtain 

2eB 2eB 
\Z{q,V,B)\^ = \Zvi{q,V,B)\^ + {n, - -^n,,re)\Z i{q.V)\\ (52) 

This expression describes quasiparticle scattering in a clean superconductor in low magnetic 
fields in a model-agnostic way, namely, it is valid regardless of the choice of the detailed model 
of quasiparticle scattering in the vortex region. Zyi{q, V, B) here describes the scattering 
off the vortex-impurity composites, which we now proceed to discuss. 



B. Impurities inside the vortex core: calculating Zyi 

As observed in the conductance ratio Z(q, V, 5), the intensity of scattering between 
parts of the Brillouin zone with the same sign of the gap grows in the magnetic field, which 
implies that the scattering potential of a vortex-impurity composite has a predominantly 
sign-preserving coherence factor. 

We now turn to a discussion of the scattering mechanisms that can enhance sign- 
preserving scattering inside the vortex cores. Table 1 shows a list of scattering potentials 
and their corresponding coherence factor effects. Weak potential scattering is immediately 
excluded. Weak scattering off magnetic impurities can also be excluded, since the change 
in the density of states of the up and down electrons cancels. This leaves two remaining 
contenders: Andreev scattering off a fluctuation in the gap function, and multiple scattering, 
which generates a t-matrix proportional to the unit matrix. 

We can, in fact, envisage both scattering mechanisms being active in the vortex core. 
Take flrst the case of a resonant scattering center. In the bulk superconductor, the effects 
of a resonant scatterer are severely modifled by the presence of the superconducting gap^. 
When the same scattering center is located inside the vortex core where the superconducting 
order parameter is depressed, we envisage that the resonant scattering will now be enhanced. 
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On the other hand, we can not rule out Andreev scattering. A scalar impurity in a d-wave 
superconductor scatters the gapless quasiparticles, giving rise to Friedel oscillations in the 
order parameter that act as Andreev scattering centersi^ii^ii^. Without a detailed model for 
the nature of the vortex scattering region, we can not say whether this type of scattering 
is enhanced by embedding the impurity inside the vortex. For example, if, as some authors 
have suggested^^, the competing pseudo-gap phase is a Wigner supersolid, then the presence 
of an impurity may lead to enhanced oscillations in the superconducting order parameter 
inside the vortex core. 

With these considerations in mind, we consider both sources of scattering as follows 



t(q, k, iun) = tA{q, k, iun) + ^^(q, k, iun) 



(53) 



where 



tA(q,k,i^„) = -Ao/A(q)(Xk+ +Xk_)Ti, 



(Andreev scattering) 



describes the Andreev scattering. Here Xk = c^; — Cy is the d-wave function with c^^y = 
cos kx,y. The resonant scattering is described by 



tif,(q>k,it^„) = 2AoSgn(u;„) //j(q)l. 



(Resonant scattering) 



Using the T-matrix approximation, we obtain for the even and odd components of Fourier 
transformed fluctuations in the local density of states due to the scattering off the super- 
conducting order parameter amplitude modulation. 



Sp' 



yj-{q,uj) = 7^ Im / Tr 



2lT 



5p"y'fici,Lu) = ^ Im / Tr 



Ck. {uj — i6) t(q, k, u; — i6) G^j^ {u — i6) 
'T3Gk_ {u — i6) t(q, k,uj — iS) {id — i6) 



(54) 
(55) 



where k-t = k ± q/2, G'k(co') = [u; — ekTs — AkTi] ^ is the Nambu Green's function for an 
electron with normal state dispersion ek and gap function Ak. We now obtain 

r even(odd) / \ £ / \ \ even(odd) / \ , _r / \ \ even(odd) / \ 

5pv (q,^) = /vi(q)A^ ^(q,a;) + /R(q)A^ ^(q,a;) 



with 



Ar"(q,^) 



...... 



2:(Ak^ + Ak_) 
{z^-ElMz^-El_) 



A--(q,a7) = ^Im 

^71" .Ik 



-i{z'^ + ek+ek_ + Ak+Ak_ 
{z^-ElW-El_) 



z=u)—i& 



(56) 
(57) 



18 



The substantially smaller odd components are: 



Ar(q,^) 



Ac 

An' 



Im / (Xk+ + Xk 



27T J I, 



-i zUk, + ek_) 



-iS 



(58) 
(59) 



where E]^ = [e^ + A^] 5 is the quasiparticle energy. The vortex contribution to the Fourier 
transformed conductance ratio fl52l) is then 



Zviiq, V, B) = nviZAiq, V, B) + Z^(q, V, 5)), 



(60) 



where 

Z^(q,V,i?) = /^(q) 
and 

Z^(q,V,5) = /K(q) 



LVo(i^) po(-V^) 



^ )Ar"(q, V) + (^ + _l_)A^^^(q, V) 



(61) 



^ ^ -)Ar"(q,v^) + (-4T + T7^)Ar(q,v^) 



>o(^) Po(-^)' 



.(62) 



IV. NUMERICAL SIMULATION 



In this section we compare the results of our phenomenological model with the experi- 
mental data by numerically computing Zyiiq, V, B) (jHUD for Andreev (1^ and resonant 
scattering. 

In these calculations we took a BCS superconductor with a d-wave gap = 
Ao/2(cosA;^. — cos ky) with Aq = 0.2t and a dispersion which has been introduced to fit 
the Fermi surface of an underdoped Ca2-xNaxCu02Cl2 sample with x = 0.12^-: 



= — 2t(cos kx + cos ky) — At' cos k^ cos ky — 2t"{cos 2kx + cos 2ky) + /i, 



where t = 1, t' = -0.227, t" = 0.168, /x = 0.486. 



A. Evaluation of Zyi 

In the absence of a microscopic model for the interior of the vortex core, we model the 
Andreev and the resonant scattering in the vortex region by constants /^(q, 2C(j„) = and 
/i^(q, iujn) = fn- Fig. 2 shows the results of calculations using these assumptions. 
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Our simple model reproduces the enhancement of sign-preserving q-vectors qi,4,5 as a 
result of Andreev and resonant scattering off vortex-impurity composites. Some care is 
required in interpreting Fig. 2, because the squared conductance ratio |Z(q, contains 
weighted contributions from both even and odd fluctuations in the density of states, with 
the weighting factor favoring odd fluctuations, especially near = 0. Both Andreev and 
resonant scattering contribute predominantly to the even fluctuations of the density of states 
(see Table 1), and give rise to the signals at ^1,4,5. In the case of resonant scattering, we 
observe an additional peak at ^3. From Table 1, we see that the Andreev and the resonant 
scattering potentials also produce a signal in the odd channel which experiences no coherence 
factor effect, contributing to all the octet q-vectors, which, however, enters the conductance 
ratio Z{ci, V) given by fl2Ul) with a substantial weighting factor. This is the origin of the 
peak at gs in Fig. 2(b). 

B. Comparison with experimental data 

The results of the calculation of the full squared conductance ratio \Z{q,V, B)]"^ are 
obtained by combining the scattering off the impurities inside the vortex core Zyi with the 
contribution from scattering off impurities outside the vortex core Zj, according to equation 
(152|1 . reproduced here: 

2eB 2eB 
\Z{q,V,B)\^ = -^ \Zvi{q,V,B)\^ + {n,--^n,ore) \Zj{q,V,B)\\ (63) 

where ricore = niTx{^/2Y is the number of impurities per vortex core. Fig. 3 displays a 
histogram of the computed field-induced change in the conductance ratio \Z{(\i,V, B)\^ — 
l-^(Qi5 ^) P cit the octet q-vectors. In these calculations, we took an equal strength of Andreev 
and resonant scattering fji = with a weak scalar scattering outside the vortex core of 
strength fi = fR = fA- In ah our calculations, we find that Andreev and resonant scattering 
are equally effective in qualitatively modelling the observations. The main effect governing 
the depression of sign-preserving wavevectors q'1,4,5 derives from the change in the impurity 
scattering potential that results from embedding the impurity inside the vortex core. 

We estimated the percentage of the impurities decorated by the vortices from the fraction 
of sample area covered by the vortices. The concentration of vortices is ny{B) = 2eB/h = 
B/^Q, where $0 = h/{2e) = 2.07 x 10~^^ weber is the superconducting magnetic flux quan- 
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(a) (b) 

FIG. 2: (Color online) Quasiparticle interference produced by the Andreev and the resonant scat- 
tering potentials, the primary candidates for producing the experimentally observed enhancement 
of sign-preserving scattering. Fig. (a) displays a density plot of the squared Fourier transformed 
conductance ratio [^^(q, predicted by ([6T|) at a bias voltage V = Ao/2 produced by pure 
Andreev scattering {/a 7^ 0, Jr = 0). Fig. (b) displays a density plot of the squared Fourier 
transformed conductance ratio \Zji{q, y)p predicted by (f62l) at a bias voltage V = Ao/2 produced 
by resonant scattering (//j 7^ 0, /a = 0). Blue lines label the positions of the sign-preserving 
q-vectors q = ^1,4,5, where both Andreev and resonant scattering is peaked. Red lines label the 
positions of the sign-reversing q-vectors q = g2,3,6,7) where both Andreev and resonant scattering 
is minimal. 

turn. The area of a vortex region is estimated as Ay = 7r(^o/2)^ with the superconducting 
coherence length ^0 = 44 A—, so that the percentage of the original impurities that are 
decorated by vortices in the presence of the magnetic field is a{B) = ny{B) Ay- Using 
these values, we obtain for the magnetic field of i? = 5 T Q;(i? = 5 T) 3.7%, and for 
i? = 11 T a{B = 11 T) 8.1%. For simplicity, we assume that a vortex core is pinned 
to a single impurity, ricore = '^i7r(^/2)^ = 1, so that the ratio of the concentrations of the 
impurities and vortices is ni/nv{B) = ricore/ Ay / {2eB /h)^ which becomes for i? = 5 T 
ml ny{B = 5T) ^ 27, and for 5 = 11 T ni/ny{B = IIT) ^ 12. 

In Fig. 3 we have modelled the scattering provided the origin of the selective enhancement 
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(a) 



(b) 



Experiment 




(c) 

FIG. 3: (Color online) Comparison between the results of the model calculations and the experi- 
mental data. Figs, (a)-(b) show the change in the squared Fourier transformed conductance ratio 
= |Z(q,F,S)|2 - \Z{ci,V,B = 0)p at q = qi-j, computed for a magnetic field of S =5 T 
(grey bars) and 11 T (red bars) at a bias voltage V = Ao/2, provided the origin of the selec- 
tive enhancement is the Andreev (Fig. (a)) or the resonant (Fig. (b)) scattering in the vortex 
core region. Here a vortex, pinned to a scalar impurity, transforms its original scattering po- 
tential with enhanced scattering at q = ^2,3.6,7 into an Andreev (Fig. (a)) or into a resonant 
(Fig. (b)) scattering potential with enhanced scattering at q = (?i,4,5 (see Table 1). Fig. (c) 
shows the experimentally observed change in the squared Fourier transformed conductance ratio 
6Z'^ = \Z{q, V, i?)p - \Z{q, V, B = 0)\'^ at q = qi-7, in a magnetic field of i? =5 T (grey bars) and 
11 T (red bars) at a bias voltage V = 4.4 meV. 



is the Andreev (Fig. (a)) or the resonant (Fig. (b)) scattering in the vortex core region. 
Both the Andreev and the resonant scattering are equally effective in qualitatively modelling 
the observations. Thus our model has qualitatively reproduced the experimentally observed 
enhancement of the sign-preserving scattering and the depression of the sign-reversing scat- 
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tering. 



V. DISCUSSION 

In this work, we have shown how scanning tunnehng spectroscopy can serve as a phase- 
sensitive probe of the superconducting order parameter. In particular, we find that the even 
and odd components of the density of states fiuctuations can be associated with a well- 
defined coherence factor. The measured Fourier transformed conductance ratio Z{q, V) = 
d//dl/(q'-l/) ^ weighted combination of these two terms, and in the limit of particle-hole 
symmetry it is dominated by the odd component of the density of states. Observation 
of coherence factor effects with scanning tunneling spectroscopy requires the presence of 
controllable scatterers. In the study by Hanaguri et al.- these controllable scatterers are 
vortices. 

Our phenomenological model of quasiparticle scattering in the presence of vortices is able 
to qualitatively reproduce the observed coherence factor effects under the assumption that 
impurity scattering centers inside the vortex cores acquire an additional Andreev or resonant 
scattering component. 

This study raises several questions for future work. In particular, can a detailed model of 
a d-wave vortex core provide a microscopic justification for the modification of the impurity 
scattering potential? One of the issues that can not be resolved from the current analysis, is 
whether the enhanced Andreev scattering originates in the core of the pure vortex, 
or from the decoration of impurities that are swallowed by the vortex core (?T-core|-Z^D/P)- 
This is an issue that may require a combination of more detailed experimental analysis and 
detailed modelling of vortex-impurity composites using the Bogoliubov de Gennes equations. 
Another open question concerns whether it is possible to discriminate between the Andreev 
and resonant scattering that appear to be equally effective in accounting for the coherence 
factor effects. 

There are several aspects to the experimental observations that lie beyond our current 
work. For example, experimentally, it is possible to spatially mask the Fourier transform 
data, spatially resolving the origin of the scattering. These masked data provide a wealth of 
new information. In particular, most of the enhancement of the sign preserving scattering 
is restricted to the vortex core region, as we might expect from our theory. However, to 
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extend our phenomenology to encompass the masked data, requires that we compute the 
fluctuations of the density of states as a function of distance from the vortex core, 

Rir, r'; ry, V) = ((5p(r - ry, l^)(5p(r' - ry, V^)), (64) 

a task which requires a microscopic model of the vortex core. 

In our theory we have used the bulk quasiparticle Green's functions to compute the 
scattering off the vortex-decorated impurities. Experiment does indeed show that the quasi- 
particle scattering off impurities inside the vortex cores is governed by the quasiparticle 
dispersion of the bulk: can this be given a more microscopic understanding? The penetra- 
tion of superconducting quasiparticles into the vortex core is a feature that does not occur 
in conventional s-wave superconductors. It is not clear at present to what extent this phe- 
nomenon can be accounted for in terms of a conservative d-wave superconductor model, or 
whether it requires a more radical interpretation. One possibility here, is that the quasipar- 
ticle fluid in both the pseudo-gap phase and inside the vortex cores is described in terms of 
a "nodal liquid"^. 

Beyond the cuprates, scanning tunneling spectroscopy in a magnetic field appears to 
provide a promising phase-sensitive probe of the symmetry of the order parameter in un- 
conventional superconductors. One opportunity that this raises, is the possibility of using 
STM in a field to probe the gap phase of the newly discovered iron-based high-temperature 
superconductors. According to one point of view^^, the iron-based pnictide superconductors 
possess an s± order parameter symmetry in which the order parameter has opposite signs 
on the hole pockets around F and the electron pockets around M. If this is, indeed, the case, 
then in a magnetic field quasiparticle scattering between parts of Fermi surface with same 
gap signs should exhibit an enhancement, while scattering between parts of Fermi surface 
with opposite gap signs will be suppressed. This is a point awaiting future theoretical and 
experimental investigation. 

We are indebted to Hide Takagi and Tetsuo Hanaguri for providing the experimental 
data. We thank Hide Takagi, Tetsuo Hanaguri, J.C Seamus Davis, Ali Yazdani, Tami 
Pereg-Barnea, Marcel Franz, Peter Hirschfeld, Zlatko Tesanovic, Eduardo Fradkin, Steven 
Kivelson, Jian-Xin Zhu, Sasha Balatsky and Lev loffe for helpful discussions. This research 
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